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Abstract 

We present investigations of the superconductor to insulator transition (SIT) of uniform a-Bi films using 
a technique sensitive to Cooper pair phase coherence. The films are perforated with a nanohoneycomb ar- 
ray of holes to form a multiply connected geometry and subjected to a perpendicular magnetic field. Film 
magnetoresistances on the superconducting side of the SIT oscillate with a period dictated by the supercon- 
ducting flux quantum and the areal hole density. The oscillations disappear close to the SIT critical point 
to leave a monotonically rising magnetoresistance that persists in the insulating phase. These observations 
indicate that the Cooper pair phase coherence length, which is infinite in the superconducting phase, col- 
lapses to a value less than the interhole spacing at this SIT. This behavior is inconsistent with the gradual 
reduction of the phase coherence length expected for a bosonic, phase fluctuation driven SIT. This result 
starkly contrasts with previous observations of oscillations persisting in the insulating phase of other films 
implying that there must be at least two distinct classes of disorder tuned SITs. 
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Superconductor to insulator quantum phase transitions (SIT) can be induced in a wide range of 
quasi two-dimensional superconducting systems, including elemental films, high Tc superconduc- 
tors, organic superconductors and superconductor graphene composites [[IHSl. Remarkably, these 
transitions occur, nearly universally, at a critical normal state resistance, -Rnc - -Rq = h/{2ey, 
and film resistances often show scaling behavior around this critical point. The most prominent 
theories that can account for these behaviors view the SIT as a Cooper pair or boson localization 
transition, rather than a Cooper pair breaking transition[|9l-[T4l. Recent experiments that probe 
films near the SIT in new ways have provided more details that support this "bosonic" picture 
of the SIT [fT5l - fT9ll . For example, scanning tunneling microscopy (STM) has revealed that spatial 
inhomogeneities develop in the order parameter on approaching the SIT suggesting that Cooper 
pairs localize into islands [fT9l l20l . High frequency transport measurements indicate that a finite 
superfiuid density persists in non-superconducting films [|21L Here, we describe experiments em- 
ploying a technique that is uniquely sensitive to the length scale that diverges at a bosonic SIT: the 
phase coherence length This technique previously revealed that Cooper pairs maintain their 
phase coherence over lOO's of nanometers through the SITs of two distinct film systems [fT5ll22l . 
For this work, we investigate a third film system [|23l - l25ll that is likely to provide the most stringent 
test of whether the bosonic SIT is generic to thin films. 

This phase sensitive technique requires patterning films with an array of nanometer-scale holes. 
The hole patterning creates a simply connected geometry that leads to Little-Parks-like (LP) 
oscillations [|26ll in the magnetoresistance (MR). The oscillations occur provided the CP phase 
coherence length {^^) exceeds the inter-hole spacing (-100 nm). For a bosonic SIT, starts at 
infinity in the superconducting phase and decreases algebraically with a power of order 1 with 
distance from the critical point [|27l l28ll . Thus, oscillations are expected to persist well into the 
insulating phase for a bosonic SIT. By contrast, in the fermionic scenario in which the order pa- 
rameter amplitude goes to zero at the SIT, collapses to zero at the critical point. We summarize 
these opposing pictures of the SIT in Fig. 1. Using this embedded hole array technique, we 
previously uncovered a CPI phase in a-Bi films deposited on anodized aluminum oxide (AAO) 
substrates[jT5l. Later we found that these substrates induced spatial film thickness variations and 
thus. Cooper pairing inhomogeneities that could give rise to a bosonic SIT ||29l . Here, we present 
results from similar experiments on a-Bi films fabricated to be uniformly thick by using atomically 
flat Si substrates. Little-Parks-like oscillations demonstrating CP phase coherence appear only on 
the superconducting side of the SIT in these uniform films. The results are consistent with 
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collapsing at this SIT in accord with a fermionic SIT. These results imply that there are at least 
two separate classes of disorder-driven SITs. 

Motivation for these investigations came from developments over the last decade that strongly 
suggest that the bosonic picture of the SIT is universal [[30l . Most convincingly, there is direct ev- 
idence that Cooper pairs (CPs) can remain intact while becoming localized and phase incoherent 
through the SIT to form a Cooper Pair Insulator (CPI) phase [|T5] - [T9l . STM experiments on super- 
conducting TiN and In oxide (InO^.) films indicated that the pairing energy gap stays finite through 
the SIT while developing spatial inhomogeneities suggestive of incipient CP islanding |fT9l l20l l3T1l . 
Other tunneling measurements on insulating InO^ films confirmed that the gap persists [fTTll . Ear- 
lier electron transport experiments on films near the SIT revealed a giant magnetoresistance (MR) 
peak larger than expected for any known single electron transport mechanism lfTSl [32] - l34l . More- 
over, strong evidence of locally phase coherent charge 2e transport in insulating InOa; films was 
obtained using a technique described below [22 J. It is now clear that this CPI phase appears in 
many film systems and its description by theories of the SIT has become pervasive. The most 
current theories of the disorder-driven SIT predict the formation of puddles of CPs that become 
phase incoherent with one another with increasing disorder[!9]-[l4l. This emergent granularity is 
predicted to appear even in uniformly disordered systems where the length scale of the disorder is 
much smaller than the correlated regions that form lfTTIl . 

On the other hand, early measurements on amorphous elemental films [|23l l24l [35] - l37]| still lend 
support to theoretical models proposing that disorder-enhanced electron electron interactions drive 
the pair density to zero at the SIT[|38l|39l. This fermionic mechanism has been used to explain the 
suppression of the superconducting transition temperature to zero with increasing disorder in ele- 
mental amorphous thin films [[38l[39l . The most significant evidence for the fermionic mechanism 
is provided by early planar tunneling data on amorphous Bi films near their SIT ll23l [351 . These 
data indicated that the CP amplitude vanishes at the transition, which creates a problem for the 
bosonic picture. However, it is still possible that the bosonic scenario is true in general if the tun- 
neling data do not reflect a vanishing CP density. The data could be explained by the appearance 
of gapless superconductivity, or a very inhomogeneous distribution of energy gaps smearing the 
density of states. In this case, superconducting pair-correlations would always exist in insulators 
close to the SIT, albeit less dramatically in some cases. 

To study the evolution of the superconducting phase coherence through the disorder-driven SIT 
of uniform a-Bi films, we quench-depo sited Inm of Sb and then subsequent layers of Bi onto a Si 
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substrate patterned with a nanohoneycomb hole array. A large area nanometer-scale hole array was 
embedded into the Si using electrochemical processes and plasma etching. A layer of AAO was 
grown into an Al film deposited on the Si using electrochemical anodization. Then a plasma etch 
was used to transfer the hole array into the Si, with the AAO as a mask. Finally, the AAO layer was 
etched away with phosphoric acid. This procedure produced an atomically flat Si substrate with 
holes in an approximately triangular array with an average hole radius of 23±12nm and spacing 
of 110±38nm (shown in the inset of Fig. 2a). The patterned Si substrate was precoated with a 
lOnm-thick underlayer of Ge and Au/Ge contact pads at room temperature. The Ge underlayer 
ensures a surface similar to a nearby reference film deposited on glass and the quench-deposited 
Sb underlayer makes the Bi wet the Si (or glass)/Ge substrate to produce uniform, amorphous 
films [|40ll4r| . Film sheet resistances, R^, were measured using standard four-point AC and DC 
techniques in the linear response regime on a (Imm)^ area of film. The sheet resistance of adjacent 
squares are similar to better than 15% at 8K, indicating good uniformity of the film thickness on 
a mm scale. Additionally, the film's normal state conductance grows linearly with thickness as 
expected for amorphous films [|29ll . 

The disorder (thickness)-tuned SIT (d-SIT) of uniform, amorphous Bi films with an embedded 
hole array (Fig. 2a) closely resembles the d-SIT of uniform, amorphous Bi films on glass (see 
e.g. Chervenak and Valles [|24l ). Thinner films are insulators as defined by dRa/dT < 0. With in- 
creasing thickness, the films become weaker insulators, developing a conductance, G = 1/R, that 
decreases as AG = aGoo lnT[? ]. For film 6, the thickest film that does not exhibit superconduct- 
ing fluctuations, a ^ I. This temperature dependence is consistent with corrections to the Drude 
conductance due to disorder enhanced electron-electron interaction and weak localization effects. 
The prefactor, a, is smaller than but of the same order as the prefactor found for films without 
holes, 1.2 ll24l . A slight increment to the film thickness produces film 7, whose R{T) turns down 
at the lowest temperatures accessed. We consider film 7, with R^ = 10.6kfi, to be very close to the 
critical point for the SIT because this downturn is a clear SC fluctuation effect ll24l . This critical 
resistance is close to, but higher than, the ^ G.Skfi value for films on glass. Subsequent evapora- 
tions produce films exhibiting sharp superconducting transitions. The transition temperature for 
SC, as measured at the resistive midpoint, grows smoothly with thickness. Thus, these films with 
holes show the same qualitative behavior as films without holes. 

Moreover, the quantitative differences in a and the critical R^ can be attributed to the influ- 
ence of the tortuous geometry induced by the hole array. The conductance of a film with holes 
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is smaller than an equivalent film without holes. The temperature dependence of the conductance 
is diminished by the same ratio. The ratio of these conductances can be calculated using finite 
element analysis of Laplace's equation for the holey geometry provided the spacing between the 
holes exceeds the electronic mean free path and the length scales that affect the quantum correc- 
tions to the conductance. For the average hole density and radii imposed by the Si substrate, this 
factor is 1.6. It is comparable to the experimentally measured ratio indicating that the microscopic 
transport properties of Bi films on holey Si and glass Il24ll or alumina substrates [f36]| are nearly 
identical. Also, it suggests that the critical normal state resistance for the SIT is dictated by the 
microscopic rather than the macroscopic properties. 

Now we turn to the detection of locally phase coherent CPs in films closest to the d-SIT. If CPs 
are present and phase coherent over a length exceeding the interhole spacing of - lOOnm, we expect 
to observe LP MR oscillations with a period corresponding to one flux quantum per hole (Hm)- 
For the hole array used to pattern these films (Fig. 2a, inset), HqHm = $o^hoies =0.22±0.02T, 
where $o = h/2e is the flux quantum and nhoies is the areal density of holes. As shown by Fig. 
2b, superconducting films (films 8 and 9) display MR oscillations on a rising background at low 
fields. The insulating films only display the monotonic rise. At high fields, the MR saturates for 
both superconducting and insulating films. Even at the lowest temperatures accessed (T=130mK) 
there are no signs of the giant MR peak characteristic of the CPI state [|T8l[33^ . [34l . 

We also measured R{T) at magnetic fields of 0, Hm/2, and Hm (Fig. 3a) to obtain the temper- 
ature dependence of the local phase coherence signal, as well as a lower temperature investigation 
of the presence of LP MR oscillations. Defining the oscillation amplitude as, 

A = i?(HM/2,T)-i?^id(T); (1) 

i?mid(T) E i (R(0,T) + i?(HM,T)) . (2) 

Fig. 3b shows that A (normalized by -Rmid) rapidly grows from zero as T drops below the mean 
field transition temperature of the superconducting films. For the film closest to the d-SIT, there is 
a hint that an oscillation due to Cooper pairing begins to emerge at the lowest temperatures. It is 
notable that this hint appears as the film resistance starts to drop. For films where no such drop in 
Rfj is observed, there is no evidence of local phase coherence. 

These results indicate that uniformly thick a-Bi films do not undergo a disorder-tuned supercon- 
ductor to CPI transition. The phase coherence length appears to collapse to zero at this SIT. While 
it is true that this behavior could indicate a first order transition at which CPs abruptly localize, 
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it seems more likely that this SIT is fermionic. Previous measurements of the tunneling density 
of states suggested that the superconducting energy gap decreases continuously toward zero on 
approaching this SIT [l23l[35l . Since the gap is proportional to the CP density, this behavior sug- 
gested that the order parameter amplitude goes to zero at the SIT. The present results support this 
interpretation. The suppression of the order parameter amplitude to zero leaves no room for phase 
coherence of the order parameter over any scale. 

To summarize, the disorder-tuned SIT of amorphous elemental films of Bi prepared with a uni- 
form film thickness involves an insulator phase that differs fundamentally from the CPI phase that 
appears for non-uniform film thickness. This interpretation implies that there are at least two dis- 
tinct classes of disorder-driven SIT; a conclusion that conflicts with recent predictions that the in- 
terplay of disorder and quantum fluctuation effects leaves A finite at the critical point [fTTl - fT3l . The 
fact that amorphous elemental films can exhibit either SIT depending on their morphology under- 
Unes the limitations of describing disorder using a single parameter like the sheet resistance. [|42ll 
Perhaps the data indicate that repulsive Coulomb interactions, which are known to reduce A with 
disorder [|38l [391 , play an essential role that varies with film structure. The resolution of this is- 
sue is important for our fundamental understanding of the SIT as well as the interpretation of 
qualitatively similar SITs, such as those of the high Tc superconductorsflJl-Sl. 
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M.D. Stewart, Jr. This work was supported by the NSF through Grants No. DMR-0605797 and 
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Figure 1: Disorder-tuned superconductor-insulator transition types Illustrations of disorder-driven SITs 
involving the emergence of inhomogeneities that lead to Cooper pair islanding (bosonic SIT) and pair break- 
ing in homogeneous films (fermionic SIT). The superconducting (SC) to insulating (INS) film transitions 
occur with increasing normal state sheet resistance, Rj^, at a critical value Rnc- Blue regions are paired, 
and pink regions are unpaired. The orange bar depicts the expected Cooper pair phase coherence length, 
which sets the size of the phase coherent regions in the bosonic insulator. The parallel stripes represent the 
phase angle of the superconducting wave function, denoted as cf). 
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Figure 2: Global and local phase coherence through the SIT a) Sheet resistance with temperature at H=OT 
for a series of uniform, amorphous Bi films quench condensed onto a Si substrate with a nanometer-scale 
array of holes (inset). The hole array has an average hole radius of 23± 12nm and spacing of 1 10±38nm. The 
film thicknesses range from 0.47 to 0.7 Inm. b) Magnetoresistance (MR) (for field applied perpendicular to 
film plane) at T=300mK for films 5-9 in (a) on a log scale. The dips spaced by Hm =0.2 IT indicate the 
presence of CPs. They appear for (superconducting) films 8 and 9 only. Inset: MR data on a linear scale for 
films 5-9 of (a) that shows the absence of an MR peakHH |3l|3l at high fields. 
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Figure 3: Temperature dependence of local phase coherence signal a) Sheet resistance with temperature 
on an Arrhenius scale for the series of films in Fig. la at three perpendicular magnetic field values: /xo-ffi =0 
(solid line), /io-f^M/2=0.105T (dashed red line), and /ioi^M=0.21T (dotted green line) (/ = H/Hm)- Inset: 
magnified view showing the drop in resistance at the lowest temperatures for film 7 at H=OT. b) Normalized 
oscillation amplitude, ^/i^mid (defined in upper inset and equations (1) and (2)), versus inverse temperature 
for films 5 and 7-9. Vertical dashed lines mark mean-field TcS for superconducting films: 0.52 and 0.78K 
for films 8 and 9 respectively. Lower inset: magnified view of the growth of oscillation amplitude at the 
lowest temperatures in film 7. 
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